Adaptive Optical Phase Estimation Using Time- Symmetric Quantum Smoothing 
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Quantum parameter estimation has many applications, from gravitational wave detection to quan- 
tum key distribution. We present the first experimental demonstration of the time-symmetric tech- 
nique of quantum smoothing. We consider both adaptive and non-adaptive quantum smoothing, 
and show that both are better than their well-known time-asymmetric counterparts (quantum fil- 
tering). For the problem of estimating a stochastically varying phase shift on a coherent beam, 
our theory predicts that adaptive quantum smoothing (the best scheme) gives an estimate with 
a mean-square error up to 2^/2 times smaller than that from non- adaptive quantum filtering (the 
standard quantum limit). The experimentally measured improvement is 2.24 ± 0.14. 

PACS numbers: 42.50.Dv, 42.50.Xa, 03.65.Ta, 03.67.-a, 06.90,+v 



Quantum parameter estimation (QPE) is the prob- 
lem of estimating an unknown classical parameter (or 
process) which plays a role in the preparation (or dy- 
namics) of a quantum system [TJ [2], and is central to 
many fields including gravitational wave inter ferometry 
[5], quantum computing [3], and quantum key distribu- 
tion [4J. The fundamental limit to the precision of the 
estimate in QPE is set by quantum mechanics [TJ [2]. 
Thus one of the key issues in QPE is the development 
of practical methodologies which allow measurements to 
approach or exceed the standard quantum limit (SQL) 
for a given measurement coupling [6] O [HJ [9j [TDJ [TTJ [12] . 
Because of its wide-ranging technological relevance, the 
prime example of QPE is estimating an optical phase 

shift [si m ma um m m h eq]. 

Apart from some theoretical papers p~9j [20] , work in 
this area of QPE has concentrated upon the problem of 
estimating a fixed, but unknown phase shift, which can 
be thought of as preparing the quantum state with an av- 
erage phase equal to this parameter. It was shown theo- 
retically [15] that for this problem adaptive homodyne 
measurements coupled with an optimal estimation fil- 
ter can yield an estimate with mean-square error smaller 
than the standard quantum limit (as set by perfect het- 
erodyne detection). This was demonstrated experimen- 
tally in Ref. [16 using very weak coherent states (for 
which the factor of improvement is at most 2). More 
recent theory and experiment have shown that interfer- 
ometric measurements with photon counting can also be 
improved using adaptive techniques [T7J [18] . 

A far richer, and in many cases more experimentally 
relevant, problem of quantum phase estimation arises 



when the phase evolves dynamically under the influence 
of an unknown classical stochastic process p~9j [20] . The 
general problem of estimating a classical process dynam- 
ically coupled to a quantum system under continuous 
measurement has recently been considered by Tsang [21] , 
who introduced three main categories of quantum estima- 
tion: prediction or filtering, smoothing, and retrodiction. 
Of those, prediction or filtering is a causal estimation 
technique that can be used in real-time applications [24] . 
Smoothing and retrodiction are acausal and so cannot be 
used in real time, but they can be used for off-line data 
processing or with a delay corresponding to the estima- 
tion time. Smoothing, in which the signal is inferred at a 
point in time based on data taken both before and after 
that time, is the only time-symmetric estimation tech- 
nique. As a consequence, it can be more precise than the 
time asymmetric techniques of filtering or retrodiction 
[20], [2T] , Such a result is very significant for quantum 
sensing applications where it is more important to have 
precise rather than real-time estimates. 

Here we present the first experimental demonstration 
of QPE using quantum smoothing. Specifically, we con- 
sider estimation of the phase of a continuous optical field, 
generalizing the theory of Ref. [19J to a more general clas- 
sical phase noise process (rather than pure diffusion), and 
to smoothing (rather than filtering). According to our 
theory, adaptive measurements and smoothing both offer 
improvements over the alternative (non-adaptive and fil- 
tering respectively). Moreover, using both together offers 
the maximum improvement, with a mean-square phase 
error smaller than the standard (non- adaptive, filtered) 
quantum limit by a factor of up to 2y/2 in theory for pure 
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phase diffusion. We verify these predicted improvements 
experimentally, for the first time in every case, and find 
a maximum improvement by a factor of 2.24 ± 0.14 over 
the SQL. 

Fig. [I] illustrates the system under consideration. The 
goal in this quantum sensing problem is to form the op- 
timal estimate Q(t) of the system phase (p(t) of a weak 
coherent state in the presence of noise in the measure- 
ment and classical noise in the system phase. The pre- 
cision of the estimate is quantified by the mean-square 
error between the estimate and the actual phase such 
thata 2 = ([e(£)-(^)] 2 ). 
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FIG. 1: Schematic diagrams showing: (a) source and local 
oscillator generation; (b) adaptive phase estimation; (c) dual 
homodyne phase estimation. Although illustrated as a single 
device, both AOMs as drawn are actually a pair of AOMs 
which shift by 110 MHz and 105 MHz in opposite directions 
to achieve a 5 MHz frequency shift. LO = local oscillator; 
RF sig = radio-frequency signal; PM = phase modulator; 
WGM = waveguide modulator; LPF = low-pass filter; MCC 
= mode-cleaning cavity; AOM = acousto-optic modulator. 

Unlike previous adaptive phase estimation experiments 
[I6j [18], we compare the phase estimate to the actual 
system phase in order to directly measure the error in 
estimation. This is achieved by deliberately imposing 
classical phase noise via an electro-optic phase modulator 
(PM), as indicated in Fig.[TJa). A titanium: sapphire laser 
operating at 860nm is used to drive the experiment. The 
arrangement of acousto-optic modulators (AOMs) shown 
in Fig.[TJa) is used to generate a pair of phase modulation 
sidebands at 5 MHz such that it is equivalent to a weak 
coherent state with a photon flux M = |a| 2 of order 10 6 
photons per second. 

The phase noise is applied using a PM driven by an 



Ornstein-Uhlenbeck (OU) [22 noise source. The phase 
variation is 



dtp(i) = -\<p(t)dt + y/RdV(t), 



(1) 



where dV is a Wiener increment and n is the inverse co- 
herence time. For this experiment k is of order 10 4 rad/s, 
so there are about 100 photons per coherence time. We 
record the imposed voltage at the monitor port of the 
high voltage amplifier (cp' in FigjlJ. The phase deviation 
arising from that voltage is calibrated via the half-wave 
voltage of the PM and used as our measure of the true 
system phase tp. 

An arbitrary quadrature of the field of interest can be 
measured with a balanced homodyne detector, in which 
the field of interest is interfered with a 1.5 mW local 
oscillator on a beamsplitter. Both outputs of the beam- 
splitter are detected and the resulting measurements are 
subtracted to form the homodyne photocurrent I(t) [23] . 
The detection efficiency (including homodyne fringe vis- 
ibility of 97%, detector quantum efficiency of 98% , and 
optical transmission of 97% ) was in excess of 89% in all 
measurements and the electronic noise floor was 11 dB 
below the shot noise of the measurements. In all cases, 
the homodyne detector is DC-locked to ensure that the 
deliberately imposed OU noise dominates the uncertainty 
in the phase. 

The adaptive phase estimation system is illustrated in 
Fig. [ljb). The output of the homodyne detector is elec- 
tronically demodulated to give I(t) which is then fed into 
the feedback filter. This yields a voltage ft \ which is 
stored for later data analysis and also fed back so as to 
imprint a phase (p oc ip' on the optical LO using a WGM. 
That is, the intermediate phase estimate cp is the phase 
of the measured quadrature. 

Because cp « (p, we can use a linearized approximation 
for the homodyne photocurrent: 



I(t)dt = 2\a\[<p(t) - (p(t)]dt + dW(t), 



(2) 



where dW(t) is Wiener noise arising from the quantum 
vacuum fluctuations. We define the instantaneous esti- 
mate 9{t) to the best estimate of cp(t) which can be made 
using only the data taken in the time interval [t, t + dt): 



9(t) := <p(t) + 

= <p(t) + 



I(t)/2\a\ 
dW(t) 1 
dt 2VM' 



(3) 
(4) 



To obtain an intermediate estimate (p with a finite 
amount of noise it is necessary to time-average the in- 
stantaneous estimate. This can be achieved by using a 
simple integrator on I(t) [19 , but for practical reasons 
we use a linear low-pass filter: 



<p(t) 



(3e ul ^ s -^I(s)/2VAfds 



(5) 
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We work in the limit where the cut-off frequency ujq « 
10 2 s _1 is much less than the feedback gain (3 > 10 5 s _1 . 
In this limit, it is easy to verify from the above equations 
that the intermediate estimate reduces to the filter used 
in Ref. [19]: 



(p{t) 



r 



(3e 



p(s-t) 



0(s)ds. 



(6) 



The intermediate estimate <p(t) of cp(t) is always a fil- 
tered estimate, because it is used in a causal feedback 
loop. The theory of Ref. [19 also used filtering to ob- 
tain the final phase estimate for ip(t). That is, the final 
estimate for time t was based only on the data obtained 
up to time t. However, it is possible to obtain a better 
phase estimate by smoothing: using the data after time 
t also. Here we can assume that data for an infinite pe- 
riod of time before and after t can be used, because the 
experimental data run time of 10 _2 s is long compared to 
the averaging time, which is < 10 _4 s. 

Let us denote by 6-(t) and 6+(t) the phase estimates 
obtained from data obtained before and after time t, re- 
spectively. Following Ref. [19] , we consider estimates that 
are weighted averages of the instantaneous estimates: 



e±(t) 



/±oo 
0(s)e Tx ^ s - t) ds. (7) 



The deviation of these from the actual phase is then 
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) dW{s). (8) 



The forwards estimate, Q-(t) corresponds to the causal 
(or filtered) estimate. A weighted average of the for- 
wards and backwards estimates can be used to construct 
the time-symmetric (or smoothed) estimate Q(t) = 
W-Q-(t) + w+Q+(i), the variance of which is: 

+2w. w+ ([e-(t) - v>(t)][e+(t) - y»(t)]>. (9) 

This can be evaluated using the definition of the OU 
process for the system phase Q. From the time symme- 
try of this process, the mean-square errors in 6+(t) and 
0-(t) are the same: 



X± 



2( X± + A) &AT 



(10) 



while the correlation between the forwards and back- 
wards estimates is 



<[e_(t)-¥>(t)][e + (t) -*>(*)]> = 



2(x-+A)( X+ + A)- 
(11) 



By symmetry, Eq. ^ expression is minimised for \- = 
X+ = X an d W- = w + = 1/2, which gives 



«(X + 2A) , X 
4(X + A) 2 



16A/" 



(12) 



Equations ( 10 ) and ( 12 ) are simplified greatly in the limit 
£ = A/ (2V kJV) <C 1, which is a good guide for our ex- 
periment whe re £ ^0.2. In this limit, the optimal value 
of x m Eq. (10) is 2^J , giving the minimum vari- 
ance g 2 _ — \JkJn /2. The relative corrections are 0(£ 2 ). 
Changing to smoothing reduces the variance by a factor 
of 2 to a 2 = 

We compare the results of the adaptive technique to 
the standard technique for phase estimation, dual ho- 
modyne detection, illustrated in Fig. [TJc). It incurs the 
same noise penalty as heterodyne detection. The dual 
homodyne data [7 + and I~ in Fig. [ljc)] can be used to 
form an instantaneous estimate, comparable to Eq. (J3|, 
via 



e B (t)=axg[I+(t)+iI-(t)]. 



(13) 



This can be shown to give an estimate which is effec- 
tively the same as Eq. Q, but with an additional noise 
penalty incurred from measuring both quadratures: M 
must be replaced by A/* s = Af/2 [15] . The mean-square 
errors for dual homodyne measurements are simply ob- 
tained by making this substitution in Eqs. JTok and (12). 



Comparing these to the adaptive results, the latter thus 
give a reduction in the variance by a factor of \/2 for both 
filtering and smoothing. The SQL is set by nonadaptive 
filtering, and equals cr 2 _ = ^k/Nq/2. 

The measured and predicted mean-square errors for 
the four different estimation techniques are shown in 
Fig. [2] The values of k and A are determined from the 
calibrated measurements of the system phase <p(i). The 
photon number, M = |a 2 |, is from the measured ampli- 
tude of the coherent state relative to the quantum noise 
limit, while x is varied in the experiment as indicated in 
Fig. [2] Each data run is 10 ms long, and error bars are the 
standard deviation of multiple data sets. We performed 
the adaptive estimates Q±(t) not by averaging 6{t) as in 
Eq. ([7]), but rather by averaging (p(t). This gives more 
stable results, and is justified since, in the limit £ <C 1, 
the optimal value of (3 in Eq. ^ is y/SxAT [19]. In the 
regime of the experiment, this (3 is much greater than 
so the extra averaging in Eq. (|6| is negligible. 

Figure [2] shows good agreement between theory and 
experiment for all the estimation techniques. It demon- 
strates that phase estimation by quantum smoothing is 
significantly better than that from quantum filtering. As 
predicted, the improvement is nearly a factor of two at 
the optimum value of x f° r both the adaptive and dual 
homodyne measurements. Figure [2] also shows the first 
experimental verification of the quantum theory of con- 
tinuous adaptive phase estimation [19]. As predicted, 
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FIG. 2: The experimental and theoretical variance a 2 of 
the four phase estimation techniques: filtered dual homodyne 
(DH) and adaptive phase (AP) in part (a); and smoothed 
DH and AP in part (b). Parameters are: kdh — 1.6218 x 
10 4 rad/s, X DH = 6.4593 x 10 4 rad/s, Mdh = 1.3235 x 10 6 
s~\ kap = 1.5868 x 10 4 rad/s, X A p = 6.1451 x 10 4 rad/s, 
Map = 1-3499 x 10 6 s" 1 . 



adaptive phase estimation outperforms dual homodyne 
measurement by a factor of approximately y/2. The the- 
ory curves here take into account the known imperfec- 
tions (detector dark noise, homodyne efficiency and op- 
tical transmission losses). However, the horizontal line 
indicating the SQL is defined (as above) in terms of 
the actual photon flux, and corresponds to what would 
be achievable by ideal dual- homodyne filtering. Note that 
adaptive measurements perform better than the SQL for 
both types of estimator. 




FIG. 3: The variance a 2 of the adaptive phase estimation for 
quantum filtering and smoothing as a function of photon num- 
ber M, compared to the relevant theoretical predictions, and 
the theoretical predictions for nonadaptive measurements. 



Figure [3] shows the optimal mean-square errors in the 
filtered and smoothed adaptive estimates for four differ- 
ent values of measured photon number, spanning an or- 
der of magnitude between M ~ 10 6 s _1 and M ~ 10 T s _1 . 
The results confirm that quantum smoothing outper- 
forms quantum filtering over a wide range of photon num- 
bers, as predicted by theory. As in Fig. [2j the SQL is set 
by ideal dual homodyne filtering, while the other the- 
ory lines are for the actual (non-ideal) experiment. As 
predicted, adaptive measurements out-perform the SQL. 

In summary, we have demonstrated experimentally 
and theoretically that estimation of the phase of an opti- 
cal field in the presence of classical noise using quantum 
smoothing is superior to the equivalent quantum filtered 
approach. We have also demonstrated experimentally 
for the first time that continuous adaptive measurements 
perform better than the standard quantum limit for both 
types of estimator. Combining quantum smoothing with 
adaptive measurements gives the maximum improvement 
over the standard (perfect non- adaptive, filtering) quan- 
tum limit. The experimental improvement of a factor of 
2.24 ± 0.14 in the mean-square error compares well with 
the theoretical maximum (using phase diffused coherent 
states) of 2y/2. These insights and techniques will be ap- 
plicable to the even more interesting case of estimation 
using non-classical states, where the improvement can be 
arbitrarily large. 
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